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ABSTRACT 



We present the work undertaken by the Evolution and 
Seismic Tools Activity (ESTA) team of the CoRoT Seis- 
mology Working Group. We have focused on two main 
tasks: Task 1 - now finished - has aimed at testing, com- 
paring and optimising seven stellar evolution codes which 
will be used to model the internal structure and evolu- 
tion of the CoRoT target stars. Task 2, still underway, 
aims at testing, comparing and optimising different seis- 
mic codes used to calculate the oscillations of models for 
different types of stars. The results already obtained are 
quite satisfactory, showing minor differences between the 
different numerical tools provided the same assumptions 
on the physical parameters are made. This work gives us 
confidence on the numerical tools that will be available 
to interpret the future CoRoT seismic data. 

Key words: Stars: internal structure, evolution, oscilla- 
tions, numerical models. 



1. INTRODUCTION 



The space experiment CoRoT will provide us with stel- 
lar oscillation data (frequencies, amplitudes, line widths) 
for stars of various masses and chemical composi- 
tions - mainly main sequence solar type stars and S 
Scuti stars - with an expected accuracy on the fre- 



quencies of ^0.1 /iHz (Ba glin et all 12002). Such an 
accuracy is needed for the determination of the po- 
sition of the external convection zone in sola r type 
stars JAudard & Provostl [1994: Monteiro eTail 120001 
l200l or the extent of the convective core of inter- 
mediate mass stars |Roxbur gh & Vorontsovt Il999t) as 
well as for the invers i on of the stellar density profile 
(Gough & Kosovichev, 1993) and of the rotation pro- 
fileor for the separation of multiplets ( Gou piTet all 
1996). However, studies in helioseismology have taught 
us that in order to probe fine details of the solar in- 
ternal structure, we need both the constraints of high 
qualit y seismic data and very accurate solar mod- 
els ( Christ ensen-Dalsgaardl 1 1 982t [Reiter et all 1 19951 
Idough et all 1 19961 iMonteiro et all I2002I) . Therefore, 
to be able to draw valuable information from the future 
CoRoT data, we have to ensure now that we will be able 
to interpret them with models having reached the required 
level of accuracy. 

Within the CoRoT Seismology Working Group, the 
ESTA group has been set up, with the aim to extensively 
test, compare and optimise the numerical tools used to 
calculate stellar models and their oscillation frequencies. 
Our goals are (i) to be able to produce theoretical seismic 
predictions by means of different numerical codes and to 
understand the possible differences between them and (ii) 
to bring stellar models at the level of accuracy required to 
interpret the future CoRoT seismic data. 

Up to now, the ESTA group has focused on specific tasks. 
The first one - Task 1 - has consisted in comparing stel- 



lar models and evolution sequences produced by seven 
participating evolution codes, presented in Section 12 In 
this context we have calculated models for several spe- 
cific study cases, presented in Section|3] corresponding to 
stars covering the range of masses, evolution stages, and 
chemical compositions expected for the bulk of CoRoT 
target stars. The main comparisons made for Task 1 
are presented and discussed in Section The second - 
Task 2 - is still in progress. It consists of testing, compar- 
ing and optimising the seismic codes by comparing the 
frequencies produced by different oscillation codes, also 
presented in Section [2] again for specific stellar cases. 
The preliminary results of Task 2 are presented in Sec- 
tion |3 The CoRoT stellar model grids are described in 
Section [6] In Sectional we present our conclusions and 
plans for future studies. 

All the results presented here are a synthesis of the 
work presented during four ESTA meetings held in 
2005, in Toulouse-France (May, 8 lh CoRoT Week 
- CW8), in Nice-France (September), in Aarhus- 
Denmark (October) and at ESA-ESTEC (Decem- 
ber, CW9). The corresponding presentations and 
posters can be consulted at the ESTA Web site at 
|http : / /www . astro . up . pt /corot /| 



2. NUMERICAL TOOLS FOR ESTA 

2.1. Stellar internal structure and evolution codes 

Seven stellar evolution codes have been used to calculate 
the models to be compared. We give here brief informa- 
tion on where to find documentation about these codes. 

(A) ASTEC - Aarhus Stellar Evolution Code: a general 
descripti on of the original code and updates can be 
found in Christensen-Dalsgaard ( 1982, 2005|). 

(C) CESAM - Code d'Evolution Stellaire Adaptatif et 
Modulair e: the original co de and updates are de- 
scribed in More]] i 19971) and lPichon & Morell ( 120051) 
and at the Website: 

|http : / / www . obs-nice . f r/ cesam/| 

(L) CLES - Code Liegeois d'Evolution Stellaire: a de- 
scription of the code, which is sti ll in an ac t ive ph ase 
of development, can be found in Scuflairel (l2005l) . 

(F) FRANEC - Pisa Evolution Code: information on the 
code and i ts updates can be found in ICariulo et al.l 
(200|) andlDe glTnnocenti & Marconilll2005l) . 

(G) GENEC - Geneva Evolution Code: the main prop- 
erties an d physical assumptions of the code are dis- 
cussed in iMevnet & Maederld2000l) . 

(S) STAROX - Roxburgh's Evolut ion Code: a. descrip- 
tion of the code can be found in lRoxburgbH2005albl) . 

(T) TGEC - Toulouse-Geneva Evolution Code: infor- 
mation on the code can be found in iRichard et alJ 
(1996) and Castr^ iEool . 



2.2. Stellar oscillation codes 

Eight stellar oscillation codes are available to us to calcu- 
late the oscillation frequencies of models. 

• ADIPLS - Aarhus Adiabatic Pulsation Package: 
available at: 

|http :/ /astro . phys . au ■ dk/~ jcd/adipack ■ n| 

• FILOU - Meudon Oscillations Code: see Suarez 
(Eool . 

• GRACO - Gra nada Oscillation Code: see 
iMovaetalJ (12004 . 

• LOSC - Liege Oscillat ions Code: see iBourv et all 
(IT?75T) : IScuflaird(l207)5l) . 

• NOC - Nice Oscillations Code: see lUnno et all 
(dH). 

• OSCROX - Roxb urgh's Oscillations Code: see 
lRoxburg"hl(l2005bl) . 

• POSC - Porto Linear Adiabatic Oscillations Code: 
See: |http : / /www . astro . up . pt/^mjm/l 

• ROMOSC - Linear N on-Adiabatic Non-Radial 
Waves: see lSuranl(ll991l) . 



3. TASK 1: INPUTS FOR MODEL COMPARISON 



3.1. Input physics for the codes 

All the models to be compared have been calculated with 
a given set of standard input physics. In this first step we 
have neglected microscopic and turbulent diffusion pro- 
cesses as well as rotation and magnetic fields. 

• Equation of State: We used the OPAL2001 
(Rogers & Navfonov, 2002) equation of state which 
is provided in a tabular form. In most codes all 
the thermodynamic quantities are obtained from 
the variables p, T, X and Z (respectively density, 
temperature, hydrogen and heavy element mass 
fraction) using the interpolation package provided 
on the OPAL Web site. But in some codes (CLES, 
STAROX) a set of thermodynamical quantities are 
interpolated in the OPAL tables (CV, P, \ p and xt 
for CLES) by a method ensuring the continuity of 
first derivatives at cell boundaries in the 4D space 
defined by the variables p, T, X and Z. The other 
thermodynamic quantities (in CLES Ti, T^ — l and 
Cp) are derived from the values of CV, P, Xp an d 
Xt using the thermodynamic relations. 

• Opacities: We us e d the OPAL95 opacity tables 
(Iglesias & Rogers, 1996) comple mented at low 
temperatures by the Alexander & Ferg uson! i 1994ft 



tables. The interpolation methods used by the codes 
may differ. In all codes the metal mixture of the 
opacity tables is fixed within a given model. 

• Nuclear reaction rates : We used the basic pp and 
CNO reaction networks up to the 17 0(p, a) 14 N re- 
action. Depending of the code the combustion of 
7 Li and 2 H is entirely followed (CLES, FRANEC) 
or these elements are assumed to be at equilibrium 
(ASTEC, CESAM, STAROX). The nuclear reaction 
rates are computed using the analytic al formulae 
provi ded by the NACRE compilation ( Angul oet all 
1999). Weak screening is assumed under Salpeter 
d!954 's formulation. 

• Convection and overshooting: We use the c l assical 
mixing length treatment of Boh m-Vitensel (1958) 
under the formulation of iHenvev et alJ i ll 9651) tak- 
ing into account the optical thickness of the con- 
vective bubble. The onset of convection is de- 
termined according to the Schwarzschild criterion 
(V ad -V rad <0) where V ad and V rad are respec- 
tively the adiabatic and radiative temperature gradi- 
ent. We adopt a mixing-length parameter q;mlt=1-6. 
In models with overshooting, the convective core 
is extended on a distance l m =a ov x mm(Hp, Rcc) 
where a ov =0.15 is the chosen overshooting parame- 
ter, Hp is the pressure scale height and i? cc is the ra- 
dius of the convective core. The core is fully mixed 
in the region corresponding to the convective and 
overshooting region; in the overshooting region the 
temperature gradient is taken to be equal to the adi- 
abatic gradient. 

• Atmosphere : Eddington's grey T(t) law is used 
for the atmosphere calculation: T=r cff [|(r+|)] 1/4 
where r is the optical depth. The level where the in- 
tegration of the hydrostatic equation starts depends 
on the codes as well as the level where the atmo- 
sphere is matched to the envelope. The radius of the 
star is taken to be the bolometric radius, i.e. the ra- 
dius at the level where the local temperature equals 
the effective temperature (r=2/3 for the Edding- 
ton's law). 

• Initial abundances of the elements and heavy ele- 



ments mixture: All models 


are calculated with the 


classical Grevesse & Noels 


(1993 


hereafter GN93) 



solar mixture of heavy elements. The mass fractions 
of hydrogen (X), helium (Y) and heavy elements 
(Z) are specified for each model. In the nuclear re- 
action network the initial abundance of each chemi- 
cal species is split between its isotopes according to 
the isotopic ratio of nuclides. 

• Physical and astrophysical constants : We used the 
reference values listed at Monteiro (2005). 



3.2. Initial parameters of the models 

We have focused on a few cases representative of the 
range of stellar masses, chemical composition and evolu- 



Table 1. Target models for Task 1 (see Fig. Q. Masses 
(M and M^ e ) are given in units of solar mass (Mq) 
while temperature is in K. "OV" indicates that overshoot 
has been included (see text). 
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Figure 1. HR diagram showing the targets for Task 1 ( see 
Table^. Red lines correspond to the PMS, black lines to 
the MS and blue lines to the post MS evolution ( CoRoT 
reference grid from CESAM). The targets are ordered in 
mass and age along the diagram, from Case 1.1 (bottom- 
right) up to Case 1.7 (top-left). 



tionary stage of CoRoT's main targets. Table ^presents 
the input parameters of the models while Figure Ogives 
their position in the HR diagram. For each model we 
examine both the evolution with time of the global pa- 
rameters (luminosity, effective temperature, mass of the 
convective core) and the internal structure parameters at 
given evolution stages. We consider several evolution 
stages: one pre-main sequence (PMS) stage with central 
temperature T c =1.9xl0 7 K, two stages close to the zero- 
age main sequence (ZAMS), two main sequence (MS) 
stages having a central hydrogen content X c of about 
half the initial value, a stage close to the core hydrogen 
exhaustion (TAMS), and a post MS subgiant stage (SG) 
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Figure 2. Evolutionary sequences and final target model 
for Cases 1.1, 1.2, 1.6 and 1.7 (all are MS). 



in which Mq b , the mass of the central region of the star 
where the hydrogen abundance is _Y<0.01, is such that 
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Figure 3. Same as Figure]2\for Case 1.3 (SG), Case 1.4 
(PMS) and Case 1.5 (TAMS). 



4. TASK 1: RESULTS OF THE COMPARISON 

4.1. Global parameters and evolutionary sequences 

For each case in Table we have compared the global 
parameters calculated by the different codes: age, ra- 
dius, luminosity, effective temperature, central tempera- 
ture and density. The results are listed in Table [2] which 
include the spread in the parameters (a.*) as defined by: 

max(cci) - min(a: l ) 

Art = 2 7 ; — 

max(Xi) + mm(arjj 



Table 2. Differences in the global parameters for all Cases (see Tabled calculated using different stellar evolutionary 
codes, [all] are the results obtained when using the global parameters from all codes; [3] when using the parameters of 
models from codes: (A, C, L) in Case 1.3, ( C, L, S) in Case 1.4; [4] is for differences between parameters of models from 
codes (A, C, L, S). 
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where the letter f3 indicates which set of codes has been 
used to calculate the model: /3=all means that values for 
all codes have been used, /3=4 is only based on values 
for codes (A, C, L, S) and (3—3 on values for codes (A, 
C, L) or (C, L, S). The subsets of codes have been selected 
because these codes have used a very similar implemen- 
tation of the physics. As can be seen in these tables when 
all models are considered in the comparisons, differences 
Ap in the global parameters may reach several percent 
and even more for the age or the mass of the convective 
core (this latter is not given in the tables). This is mainly 
explained by the fact that the reference set of physics has 
not been fully implemented in all codes. For the set of 
codes that have followed more closely the physics speci- 
fied, the differences are significantly smaller as illustrated 
by A4 and A3. The corresponding evolutionary tracks 
in the HR diagram plotted in Figures |2] and [5] are fairly 
consistent as expected from the comparison of the global 
parameters of the target models. 

Several iterations have been made in the comparison pro- 
cess. This has allowed us to find some errors or inaccu- 
racies in the evo lution codes or in the model calcu lation 
proce s s (se e Monteiro e{ al||2 005tlL ebreton & Monteiro, 
| 2005t IChristensen-Dalseaardl 120051 tMontalba n et all 
2005a b, and references therein). As a result, the relative 
differences between the models have decreased. How- 
ever, further work has to be done for the codes which still 
do not follow exac tly t he physics specifications decided 
for Task 1 (Section l3~H and for which the greatest differ- 
ences with the other codes are found. 



4.2. Internal structure 

We have compared the variations of the hydrogen abun- 
dance X and of the square of the sound speed c 2 inside 
models calculated by different codes. We have only used 
models whose global parameters are very similar in or- 
der to ensure that the differences are mainly determined 
by how each code calculates the evolution and the struc- 
ture of the model and not by significant differences in the 
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Figure 4. Plot of the relative sound speed differences and 
hydrogen abundance differences at fixed relative mass be- 
tween pairs of models corresponding to Case 1.5 (with 
overshooting). 



physics. Some model differences in X and c 2 are shown 
in Figures 0JJ5] while Table [3] provides the spread of the 
model differences in the interior. The output from CE- 
SAM has been used as the reference for these differences. 
In this comparison we have excluded the outermost sur- 
face region (to>0.95 M) where we have found substan- 
tial differences between models; understanding the origin 
of these differences requires further work. 




Table 3. Upper limits for the model differences in sound 
speed squared ( c 2 j and in hydrogen abundance (X) for 
the interior structure (m/M <0.95) of the models being 
compared. The CESAM model is used as the reference for 
all cases while the last column indicates the codes used 
for the comparison. These values are a summary of the 
interior differences shown in Figures^^ 
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8%) occur in the region of the convective core. Models 
differ in the treatment of the overshooting and of the con- 
vective boundaries and differences are spread in a mass 
range between ^0.05 and 0.20 m/M because the extent 
of the convective core decreases with evolution. 

The sound speed profile comparison also shows differ- 
ences in the central regions as well as some differences 
in more external regions. These reflect differences in the 
treatment of the physics, in particular the equation of state 
and its implementation in the codes. The representation 
of the network of reactions is also the source of some dis- 
crepancies. Figure |5] illustrates the need to work further 
on how the different species are treated in the reactions 
network and a more consistent use of the equation of state 
tables, although the differences for these Cases are much 
smaller (below 0.1% for X and below 1% for sound- 
speed squared). A factor in explaining the differences in 
the pre-main sequence phase is that different models are 
started with different initial conditions. 



One of the key aspects that dominates the model differ- 
ences is the treatment of the borders of convective re- 
gions and how these change in time. This is clear in 
Case 1.5 (Figure |4j, corresponding to the 2 M Q model 
on the TAMS with overshooting, large differences (up to 



These are some of the items that must be further stud- 
ied in order to overcome part of the remai ning model 
differenc es. Also, as demonstrated by Montalba n et alJ 
(2OO5a0) from the comparison of the CLES and CESAM 
codes, differences can be reduced to a very small level 
(generally less than 0.5%) when similar prescriptions are 
adopted for the input physics and their implementation. 
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Figure 6. Plot of the frequency differences, between mod- 
els produced by different codes, for Case 1.2. The full line 
is for £=0, dotted line for £=l, dashed line for £=2 and 
dot-dashed line for £=3. Also shown are the frequency 
differences when the scaling due to the stellar radius (R) 
is removed. The remaining differences in the frequencies 
are mainly due to near surface effects. 

4.3. Seismic properties 

For the comparison of the seismic properties of the mod- 
els we have used one seismic code POSC in order to cal- 
culate the frequencies of oscillations of all models for 
Case 1.2. We consider modes with degree £=0,1,2,3 
and order n=5— 20. In Figure[6]the frequency differences 
are shown for these modes. After removing the scaling 
due to different stellar radius, the differences of the fre- 
quencies are dominated, as expected, by the near surface 
differences of the models. The component due to dif- 
ferences in the deep interior, as given by the differences 
between curves for different mode degrees, is small, be- 
ing below the limit of the expected accuracy of the fre- 
quencies to be provided by CoRoT. For the same models 
in case 1 .2, we also compared the frequency separations, 
i.e. the large frequency separation A^=v n+ i^—i/ nt g be- 
tween the frequencies (y n i) for modes of degree £=0 
and consecutive mode order n and the small frequency 
separation 8i : i+2=Vn,i—v n -i,i+2 between the frequen- 
cies (y n i) for modes of degree £=0 and £=2 with mode 
orders n and n—1 respectively (see Figure^. The calcu- 
lated separations are very similar, in accordance with the 
similarity of the global parameters of the models. 
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Figure 7. Plot of the large ( Aq) frequency separation for 
Case 1.2. The small ( 602) frequency separations are also 
shown. 



5. TASK 2: FREQUENCY COMPARISON 



In order to test the seismic codes being used we have 
started by comparing the frequencies of p-modes calcu- 
lated using the same model (at exactly the same fixed 
mesh). For this exercise a 1.2 Mq stellar model was 
used. The spectra considered in the comparison include 
frequencies of modes with degree £=0, 1,2 and 3. In 
all comparisons we use LOSC data as the reference. 
Some codes have used Richardson extrapolation and oth- 
ers have not. 

Figure[8]shows the relative frequency differences for seis- 
mic data calculated with different codes (and different nu- 
merical procedures). Differences for £—0 modes present, 
in general, larger values than for 1^0. The figure also 
confirms that the use of Richardson extrapolation reduces 
the discrepancies between codes. 

Figure [9] shows relative differences for the large sepa- 
ration (top) and small separation (bottom) for the same 
codes and spherical degrees as in Figure [8] Codes us- 
ing Richardson extrapolation present very similar results, 
while £—0 have larger differences than those coming 
from 1=1. The relative differences of the small separa- 
tion 5o2 are one order of magnitude larger than 5\z- 

A summary of the results is listed in Table |4] The major 
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Figure 8. Relative frequency differences as a function of 
the frequency for modes with £=0 (dots) and £=1 (lines). 
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Figure 9. Top: relative differences of the large sepa- 
ration as a function of the frequency for £=0 (dots) and 
£—1 (lines) modes. Bottom: relative differences of the 
small separation 6q2 (lines) and <5i3 (dots) as a function 
of the frequency for modes. 



source of the differences is the use, or not, of Richard- 
son extrapolation, but different codes react differently 
to the mesh as some are second order schemes while 
other used fourth order integration schemes (see also 
iRoxburghl l2005bl) . A further step on the analysis on the 
sensitivity of the codes to the mesh of the models and 
the treatment of th e boundary conditions is in progress 
dMova et al.l [20051) . 



Table 4. Differences between frequencies and differ- 
ences of the frequency separations calculated using the 
oscillation codes listed in Section^\for a stellar model of 
1.2 M Q at the ZAMS. 
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6. REFERENCE GRIDS 



In parallel to Task 1 and Task 2, grids of models have 
been especially calculated with the CESAM and CLES 
codes for masses in the range 0.8—10 Mq and chemical 
compositions [Fe/H]=0.0 and —0.10. These reference 
grids have already been used to locate C0R0T potential 
targets in the HR diagram and study some 5 Scuti candi- 
dates for C0R0T ( Poretti et al., 2005). In addition oscilla- 
tion frequencies for stellar models from the CESAM code 
have also been calculated either with POSC or ADIPLS 
codes. This material can be found on the ESTA Web site 
at |http : / / www . astro . up ■ pt / cor ot /model s/| 



7. CONCLUSIONS AND PERSPECTIVES 



The comparisons undertaken under Task 1 have shown 
that stellar models calculated by seven different codes are 
to first order consistent. Differences in the global param- 
eters and in the internal hydrogen abundance and sound- 
speed profiles are generally in the range 1-10 percents 
except for the age and the extent of the convective core 
that may differ by 30% or more. 

For codes that have followed precisely the reference set 
of physics defined for the exercise, the agreement is much 
better with differences of the order of 1 percent or less, 
which gives us confidence for future modelling of C0R0T 
targets by these codes. Nevertheless, even if the refer- 
ence set of physics is fixed, differences in the models re- 
main due to different implementations of physics in the 
codes, different ways of interpolating in the tables, dif- 
ferent choices of input parameters (isotopic ratios, atomic 
masses), etc. In a detailed comparison o f the CLES and 
CESAM codes, Montalba net al.l (|2005a b) have tried to 
identify these remaining inconsistencies showing that it 
should be possible to limit the differences to the order of 
0.5% or less. 

Concerning the global parameters, the source of the 
sometimes still large differences in age have to be exam- 
ined which will imply the adoption of a common defini- 
tion for the age zero. Concerning the seismic parameters, 
the differences in the total radius and in the structure of 
the surface layers are responsible for the main differences 



in the frequencies, thus requiring a more detailed compar- 
ison of these zones. 

Task 2 must be further developed before we can assess if 
the seismic codes are providing precise frequencies. Fur- 
ther work on the mesh of the model, the numeric scheme 
used to solve the oscillations equations and the boundary 
conditions implemented is in progress. 

A detailed evaluation of the numeric precision and of the 
implementation of physics in the codes will continue to 
be necessary as the physics is improved and the methods 
and assumptions in the calculation become more com- 
plex. An example is the inclusion of diffusion to calculate 
the evolution - how this is implemented requires a de- 
tailed study. Consequently, the work initiated by ESTA 
will proceed after CoRoT flies. As long as we need to 
calculate new stellar models with better physics, the need 
for a regular consistency analysis of the codes remains. 
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